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Abstract 

In this paper we analyze the behavior of quantum random walks. In particular we present 
several new results for the absorption probabilities in systems with both one and two absorbing 
walls for the one-dimensional case. We compute these probabilities both by employing generating 
functions and by use of an eigenfunction approach. The generating function method is used to 
determine some simple properties of the walks we consider, but appears to have limitations. The 
eigenfunction approach works by relating the problem of absorption to a unitary problem that 
has identical dynamics inside a certain domain, and can be used to compute several additional 
interesting properties, such as the time dependence of absorption. The eigenfunction method 
has the distinct advantage that it can be extended to arbitrary dimensionality. We outline the 
solution of the absorption probability problem of a (d-l)-dimensional wall in a d-dimensional 
space. 



1 Introduction 



Several recent papers have studied the properties of quantum walks, which are quantum computa- 
tional variants of discrete-time random walks & [|, |8[ ^, 11, 12, 13, 14, 15, [II], 21]. The behavior 
of quantum walks differs from that of ordinary random walks in several striking ways, due to the 
fact that quantum walks exhibit interference patterns whereas ordinary random walks do not. For 
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instance, mixing times, hitting times, and exit probabilities of quantum walks can differ signifi- 
cantly from analogously denned random walks g, ||, [|, 15]. One-dimensional quantum walks are 



also relevant to quantum chaos fl2Q| . (Continuous-time variants of quantum walks have also been 
proposed and exhibit significant differences from classical continuous random walks ||, 0- However, 
we will only discuss the discrete-time case in this paper.) 

Ordinary random walks have had many applications in computer science, particularly as algo- 
rithmic tools. Examples include randomized algorithms for graph connectivity, 2SAT, and approxi- 
mating the permanent (see, for instance, Lovasz [10| for several examples of algorithmic applications 
of random walks). Quantum walks have the potential to offer new tools for the design of quantum 
algorithms, which is one of the primary motivations for studying their behavior. For example, 
given that quantum walks on certain simple structures mix faster or have faster hitting times than 
random walks, there is the potential that they will admit quantum speed-ups for algorithms based 
on random walks on more complicated structures. 

In this paper we investigate the behavior of a class of one-dimensional quantum walks that are 
very simple generalizations of the quantum walks introduced in Refs. 12, |3). 



One-dimensional quantum walks are not likely to be directly applicable to algorithm design. 
Although there has been some work on quantum walks on general graphs [0], many questions 
about quantum walks on general graphs appear to be quite difficult to answer at the present time. 
Thus it is important to extend techniques developed to analyze one-dimensional quantum walks to 
quantum walks on general graphs. To this end, we also investigate the behavior of quantum walks 
in d dimensions. 

Ref. H demonstrates that one-dimensional quantum walks differ qualitatively from classical 
random walks. For example, this walk spreads with time as t instead of \[t. Moreover, if the 
walk is evolved in a system with one absorbing boundary, the probability of eventual absorption 
by the wall is less than unity (in contrast to a classical unbiased one-dimensional random walk, 
where this probability is unity). In this paper we generalize and extend the results of Ref. Q and 
calculate the absorption probabilities of one-dimensional quantum walks in systems with one and 
two absorbing walls. We use a combinatorial method as well as an eigenfunction expansion method. 
The combinatorial method is similar to that used in Ref. || , and the eigenfunction method has been 
used by others to address periodic systems ||, [U|, [l6|, |2[ as well as systems with potential 
steps [13, 14]. Naturally the two methods agree perfectly in all cases in which results have been 
obtained using both methods. We extend the eigenfunction method to general dimensionality. 

Our results may be summarized as follows. First, for one-barrier systems, we obtain exact 
expressions for the probability of absorption by the barrier, as a function of the initial distance 
to the barrier. (Complementing these formulas gives the probability of escape to infinity.) These 
expressions involve integrals of different forms coming from our two methods of analysis. Both forms 
allow asymptotic analysis of the absorption probabilities; in particular, we compute the limiting 
probabilities when the initial distance to the barrier is large. We do this both for Hadamard walks 
and for walks based on more general unitary transformations. Next, for the two-barrier Hadamard 
system, we compute the long-time limit of the probability of absorption by each barrier when 
the walker starts off very far from one barrier but an arbitrary distance from the other barrier. 
Again, the expressions involve integrals whose asymptotic limits are easily analyzed, so that we 
can compute the limiting probabilities when the initial distance to both barriers is large. We then 
outline how the eigenfunction method can be used to analyze the behavior in small systems. Then 
we use the eigenfunction method to analyze the time dependence of the absorption in the limit of 
long times for walks with both one and two walls. We find that the approach to the asymptotic 
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limit is much slower when there are two walls, where the probability remaining to be absorbed at 
time t decays as than when the system has one wall, where the probability remaining to be 

absorbed at time t is proportional to 1/i 2 . Finally, we study d-dimensional walks, and show that 
the region over which most of the probability is distributed by time t has volume proportional to 
t d . We indicate how to solve the problem of the absorption of a (d — l)-dimensional wall without 
giving explicit results. 

The paper is organized as follows. In Section [| we give definitions for one-dimensional quantum 
walks and for the specific processes based on quantum walks that we consider. Section |3] presents 
our results based on the combinatorial approach and Section || presents our results based on the 
eigenfunction method. Finally, in Section |B| we analyze generalizations of one-dimensional quantum 
walks to higher dimensions. 

2 Definitions 

2.1 One-dimensional quantum walks 

For any finite or countable set S we may denote by 7~i(S) the Hilbert space of all square-summable 
functions from S to the complex numbers C, along with the usual inner product. Using the Dirac 
notation, the space H(S) has a standard basis {\s} : s G S}, which is orthonormal. One-dimensional 
quantum walks are discrete-time quantum processes on the space W(Z x {L,R}). The standard 
basis for this space therefore consists of elements of the form \n,d), where n £ Z is the location and 
d 6 {L, R} is the direction component of such an element. 

Given an arbitrary unitary operator U on 7i({L, R}), define a unitary operator W v acting on 
H(Z x {L,R}) as follows. For each standard basis state \n,d) we have 

Wu\n,d) = (L\U\d) \n-l,L) + (R\U\d) \n+l,R) , 

and we extend W v to all of 7i(Z x {L, R}) by linearity. Alternately we may define W a = T(I ® U), 
where T is defined by 

T \n, L) = \n — 1, L) , T \n, R) = \n + 1, R) , 

and we identify 7i(Z x {L,R}) with the tensor product space (g) TC({L, R}) in the natural 

way. We use the term one- dimensional quantum walk to refer generally to any process involving 
the iteration of W n , since such processes are reminiscent of a particle doing a random walk on a 
one-dimensional lattice. 

We will also consider quantum walks on higher dimensional lattices in the final section of this 
paper — definitions for this type of walk appear in that section. 

2.2 Absorbing boundaries and exit probabilities 

In this paper we will be interested in the situation in which our system is initialized to some state 
and we alternately apply the operator W v and perform some measurement of the system. The type 
of measurements we focus on are as follows. For each n S Z, consider the projections IIy es and II™ 
defined as 

n" s = \n, L) (n, L\ + \n, R) (n, R\ , Il£ = I - II" . 
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These projections describe a projective measurement that corresponds to asking the question "is 
the particle at location n?". Given a system in state \ip), the answer is "yes" with probability 
||IIy CS \ ip) || 2 , in which case the state of the system becomes n™ es \ tp) (renormalized) and the an- 
swer is "no" with probability ||II™ \tp) || 2 , in which case the state of the system becomes H™ 
(renormalized). 

The first type of process we consider is the one-boundary quantum walk, which is as follows. 
The system being considered is initialized to some state |0) (a \L) + J3 \R)), which corresponds to 
a particle at location and having direction component in state (a \L) + f3\R)). Fix an integer 
M > 0, which will be the location of our absorbing boundary. For given U, we alternately apply Wu 
and the measurement described by {n.y' es , n* 7 }, which gives result "yes" if the particle has reached 
location M and "no" otherwise. The process is repeated until the result "yes" is obtained. The 
probability that the result "yes" is obtained is the exit probability for this walk. For fixed U and 
given M, a, and /3, we will denote this exit probability by r M (a,(3). We also write p M = r M (l,0) 
and q M = r M (0, 1) for short, i.e., p M is the exit probability for starting in state |0, L) and q M is the 
exit probability for starting in state |0, R). 

We also consider two-boundary quantum walks. In this case, the particle is initialized in some 
state |0) (a \L) + j3 \R)) and we alternately apply W a and the (commuting) measurements described 
by {II'Ml ,n-*' lL } and {ltff e f ,I^f R } for < M L ,M R . The quantities of interest in this case are 
the probability of exiting from the left (i.e., measuring the particle at location —Ml) and the 
probability of exiting from the right (measuring the particle at location Mr). Again, in this case 
the boundaries are absorbing, since the process is terminated when either measurement gives result 
"yes". 

2.3 Statement of results 

A summary of our results is as follows. 

First, we give a complete description of the exit probabilities for one-boundary quantum walks 
in one-dimension for arbitrary unitary U, any starting state of the form |0) (a\L) + (3\R)), and 
any boundary location M > 0. This includes integrals for computing exact exit probabilities, 
closed form solutions for the exit probabilities in the limit of large M (one of which proves a recent 
conjecture of Yamasaki, Kobayashi and Imai [^]), and several other results concerning the behavior 
of these probabilities. It is proved that it is enough to consider only real unitary U for the purposes 
of analyzing such walks. 

Next, for two-boundary Hadamard quantum walks in one dimension, we present integrals for 
the exact exit probabilities when the particle starts out an arbitrary distance from one wall and 
an asymptotically large distance from the other. These integrals are evaluated to yield closed 
form solutions for the exit probabilities when the distance from the first wall is small and when 
the particle starts out very far from both walls. We then calculate the time-dependence of the 
absorption probability at long times for both one- and two-wall walks. 

Finally, for d-dimensional walks we generalize some of the results of Ref. |J on one-dimensional 
walks, including a derivation of the asymptotic form for the amplitudes associated with the walk 
and a demonstration that the amplitude spreads nearly uniformly. For (f-dimensional quantum 
walks with a (d — l)-dimensional barrier an integral for the absorption probability is derived. 

3 Combinatorial analysis 
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3.1 Generating functions 

As described in Section [2], we consider the case where our system is initialized in some state 

|0) (a\L) + p\R)) 

and we alternately apply the operator W v and measurement given by {n^L,, 11^} for some M > 0. 
This process is repeated until the measurement gives result "yes", at which time the process is 
terminated. 

We begin with two special cases: the first is the case that the starting state is |0,L) and the 
absorbing boundary is at M = 1, and the second is the case that the starting state is |0, R) and 
the boundary is at M = 1. We will define generating functions for these cases that are used to 
determine exit probabilities for all starting states and boundary positions. For given unitary U 
define two generating functions / and g as follows: 

oo 

fu(z) = < M > R \ WviKoWvf' 1 |0, L) z\ 
t=i 

oo 

g v (z) = J2 < M > R \ WvtP&W,,)*- 1 |0, R) z\ 
t=l 

The coefficient of z in fu( z ) is therefore the (non-normalized) amplitude with which the system is 
in state \M,R) after t time steps, assuming the system starts in state |0,L), and gu{z) is similar 
except we start in state |0, R). We will simply write / and g to denote f u and g u when U is 
understood. Thus, for example, the probability that a particle starting in state |0, L) is eventually 
observed at location is 

pi=EI^]/(*)l 2 > 

t>0 

where denotes the coefficient of z l in f(z), and similarly the probability that a particle 

starting in state |0,i?) is eventually observed at location is 

t>0 

The reason that the generating functions / and g are useful for analyzing exit probabilities for all 
boundary positions is as follows. For given M > 2 consider a generating function defined similarly 
to /, except for the boundary at location M rather than location 1. Then this generating function 
is simply f(z)(g(z)) , which follows from the fact that to get from location to location M, the 
particle needs to effectively move right M times, and for each move after the first, the direction 
component is R. Similarly, the generating function corresponding to starting in state |0, R) is 
simply (g(z) ) M . 

Given arbitrary generating functions u and v, their Hadamard product is u © v, defined by 

(u®v)(z)=J2(lz t Hz))([z t ]v(t))z t . 

t>0 

Thus, pi = (/ /)(1) and q\ = (g In general we have 

(u „)(1) = i- j 2 J u(e* e )v(e~ w )de, (1) 
provided J2t>o([ zt ] u ( z ))(l zt ] v (t)) converges. (This follows from results in Section 4.6 of |jl8|1 .) 
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3.2 Hadamard Walk 



The most common choice for U in recent papers on one-dimensional quantum walks has been (or 
is equivalent to) the following: 

U \L) = ^(\L) + \R)), U \R) = -L(\L) - \R)), (2) 

i.e., U is the Hadamard transform where we identify \L) = |0) and \R) = |1). The resulting walk 
has been called the Hadamard walk. It turns out that the general behavior of this walk is not 
specific to the Hadamard transform, as we will show shortly. (Nayak and Vishwanath [16] have 
also claimed results concerning the generality of the Hadamard walk.) However, it is helpful to first 
consider the Hadamard transform because it is simple and we can reduce the behavior of general 
quantum walks to the Hadamard walk. 

It can be shown that for U as in Eq. (pf) we have 



/(*) = — 



V2z 

i-z 2 - VTTz^ ( , R 

g(z) = = f(z) - V2z. 

We will not argue this here, since later we will derive generating functions for arbitrary U that give 
these generating functions in the case of the Hadamard transform. 
Define F(0) = f(e id ) and G{9) = g{e ie ). Then by Eq. (|l|) 

Pi = ± / \F(8)\ 2 d9 = - 

(■2tt 



Qi = [ 7T \G(9)\ 2 d9 = -, 
2ir J vr 



as proved in Ref. ||]. Thus, in this case a particle starting at location 1 has a 1 — 2/7T ~ 0.3634 
probability of "escaping" the absorbing boundary at location 0, which contrasts with the classical 
unbiased random walk, for which the probability of escape is 0. 

Suppose now that the boundary is at location M for any M > 1. Then from the discussion in 
the previous subsection we conclude that 

Pm = ^ \F(9)\ 2 \G(9)\ 2M - 2 d9 (3) 



q M = ±- j 2 J \G{9)\ 2M d9. 



These expressions make it very easy to calculate the exit probabilities in the limit for large M. 
Since \G(9)\ < 1 for all 9 G [0, 2vr], with strict inequality for 9 £ (vr/4, 3vr/4) U (5vr/4, 7vr/4), we have 

= kf> w?dM+ lf^ Fm2M 

4 4 

2 1 

= ps .1366 

vr 2 
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and 

•■■'It, 

1m = lim — / 



lim q M = lim — I \G(6)\ 2M dO 



1 /"4 1 — 

= — d6 + — d9 

4 4 

1 

2' 

Now let us calculate the exit probabilities in the limit for large M for arbitrary directional 
component in the starting state. Assume the particle starts in state 

|0> {a\L)+0\R)). 

Recall that we denote the exit probability for starting in this state by r M (a,(3). The generating 
function for absorbed paths is now 

af(z)g(z) M - l +(3g(z) M = (af(z) + [3g{z))g{z) M ' 1 . 

The exit probability is therefore 

1 f 27r 

r M (a,/3) = — \aF{9)+f3G{e)\ 2 \G{6)\ 2M - 2 d6 
^ Jo 

~ 77- / \aF(6)+f3G{d)\ 2 d6 
2tt J t 



where T = (-f , f ) U (2f , ^f). It will be helpful to note that 



1 _ 1 1 

7T 7T 2 ' 



which follows from the fact that g(z) = f(z) — y/2z. Therefore 



This probability is maximized when a = sin(7r/8)e 1 ^ and (3 = — cos(7r/8)e^, (for arbitrary real 4>) 
giving 

2 \/2 

lim r M (a,(3) = —= H .5753 

and is minimized when a = sin(37r/8)e 1 ^ and (3 = cos(37r/8)e"^, giving 

lim r M (a,(3) = — \= + 1 + — » .0614. 

M^oo a/2 7T 
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3.3 Asymptotics of Exit Probabilities 

It is clear from the integral representation of r M (a, (3) and the estimate \G\ < 1 that these probabil- 
ities decrease with M, to the limits computed in the last section. In this subsection we investigate 
how quickly the limit = lim M _^oo r M (a, (3) is approached. We shall prove that 

r M (a,(3)=r oo + 0(M~ 2 ) 

and give an asymptotic series for the remainder. The principal technique is Watson's lemma Q, 
which in this case reduces to successive integration by parts. 

Expanding the integral representation of the last section, we get 

r M (a,P) = !S£ \F\ 2 \G\ 2M ~ 2 dO 
27r Jo 

+ Ifi! [ 27T \G\ 2M d0 

+ [ 2 \(f/g)\g\ 2M d9. 

2?r Jo 

(The last term has no contribution from 9 (a/3) because of symmetry.) 
The coefficient of |ct| 2 is 

p M = JL j 2 J \ F \ 2 \G\ 2M ~ 2 de = i- fj \F/G\ 2 (|G| 2 ) M de. 

As a function of 9, the integrand is analytic with the possible exception of branch points occurring 
when 8 = ±7r/4, ±37r/4. Explicitly, when 7r/4 < 9 < it/ 2 we have 

2 _ (2 sine - V-2cos 29) 2 
\G\ - 

and 

\F/G\ 2 = 1 - 2 cos 26 + 2 sin 6^-2 cos 20", 

as can be seen from the formula \a + 6| 2 = \a\ 2 + |6| 2 + 23f?(a6). 

Making the substitution e _n = |G(#)| 2 and integrating by parts v times gives 

r /2 f log(3+2v / 2) " d k UO) 

Iu ■= \F\ 2 \G\ 2M ~ 2 d6 = jf ^u)e^ M du ~ £ ^ + 0(l/M^ 2 ). 

We now determine the power series for the multiplier <\> around u = 0, as follows. Observe that 
when 6 = tt/4 + t, we have 

(2 sin (9- V-2cos26>) 2 , „ /= a2 ^ , 
- ^ £- = 1 - 2V2s + 4s 2 - 2V2s 3 + • ■ • 

with s = t 1 / 2 , whereas 



By the implicit function theorem, there is a function r, analytic around 0, for which s = r(u). The 
coefficients of r can be computed term by term; we have for example 

, x y/2 y/2 3 

r[u) = — u H u + • • • . 

w 4 96 

Since r'(u)du = (l/2)t~ 1 / 2 dt, we must have 

cj>{u) = 2r{u)r'{u)\F/G\ 2 8=r{u) . 

Coefficients of this may be found using r and the expansion of \F/G\ 2 as a power series in s. 
Explicitly we have 

xl n 1 1 2 
The symmetry and periodicity of the integrand imply that 

'2 i\ i r/ 2 



From this and the above we get the asymptotic series 

(2 _ 1\ 1 1 2 4 79 

Pm ~ V ~ 2/ + 2^M2 + ^M3 + ^4 + ^5 + + ' ' ' W 

Aficionados of the "law of small numbers" will appreciate that one must develop this series to order 
6 to obtain a coefficient that does not fit the initial pattern. 
Applying similar reasoning, the coefficient of \(5\ 2 is 

1 1 1 19 

9M ~2 + 2^ + 27M^ + 8^ + ''' ' (6) 

and the coefficient of 25ft (a/3) is 

(1 _ 1\_ 1 _ 3 _ _2_ _ 15 

3.4 Other Transformations 

In this section we argue that the exit probabilities of the Hadamard walk are not really specific to 
the Hadamard transform. The argument can be generalized to other properties of the Hadamard 
walk. More generally, we show that it suffices to analyze unitary transformations U with only real 
entries. 

Suppose instead of using the Hadamard transform we let U be the general transformation 
defined by: 

U\L)=a\L) + b\R), U \R) = c\L) + d\R) . (8) 

We will consider generating functions f v (z) and Qu{ z ) defined in Section [T^ for this general trans- 
formation U. It is easy to see that these generating functions must satisfy 

f v (z) = bz + azfu(z)g v (z) 
g v (z) = dz + czf u (z)g u (z). 
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Solving these equations for f v (z) and gu(z) and taking the solutions that make sense for small 
powers of z gives 



fu(z) 



9u(z) 



I -{ad- bc)z 2 - y/\ - 2(ad + bc)z 2 + (ad - bc) 2 z 4 



2cz 



1 + (ad - bc)z 2 - y/1 - 2(ad + bc)z 2 + (ad - bc) 2 z 4 

2az 



The first few terms of these functions are as follows: 

fu(z) = bz + abdz 3 + abd(ad + bc)z 5 + abd(a 2 d 2 + 3abcd + b 2 c 2 )z 7 + 
abd(a 3 d 3 + 6ab 2 c 2 d + 6a 2 bcd 2 + b 3 d 3 )z 9 + ■■■ 

9u ( z ) = dz + bcdz 3 + bcd(ad + bc)z 5 + bcd(a 2 d 2 + 3abcd + b 2 c 2 )z 7 + 
bcd(a 3 d 3 + 6ab 2 c 2 d + QaHcd 2 + b 3 d 3 )z 9 + ■■■ 

Letting X = ad and Y = be, we see that 



fu(z) = 



9u(z) 



1-(X- Y)z 2 - y/1 - 2(X + Y)z 2 + (X - Y) 2 z 4 



2cz 



1 + (X- Y)z 2 - y/1 - 2(X + Y)z 2 + (X - Y) 2 z A 



2az 



Next, we can simplify matters by taking into account that U is unitary An arbitrary 2x2 
unitary matrix can be written 

e^+^^/p e^-t+^y/T^p 
eW-My/T^p -e^-^y/p 



(9) 



where r), 4>, if) are real and < p < 1. The global phase of r\ will not affect the behavior of the walk 
in any way, so for simplicity we may set rj = without loss of generality. We will write U p ^^ to 
denote this transformation in order to stress the dependence on p, <f>, and ip. This leaves us with 

X = -p 
Y = I -p. 

Since X — Y = — 1 and —(X + Y) = 2p — 1, we may write 

1 + Z 2 - y/l + 2(2p-l)z 2 + Z 4 def 



9u a ,s.Az) 



2e i (~<f> + ^ y/T^pz 



1-z 2 - y/1 + 2(2p - 1) z 2 + z 4 
2e^^y/pz 



def 



- 9p,<p,i>{z)- 



Thus, 



\F P ,4,^(6)\ 2 
\Gp, ( j >)1 p(6)\ 2 



1 



4(1 - p) 



1 + e 2ie - yjl + 2(2/> - l)e 2id + , 



1 

Tp 



l-e 



2w 



1 + 2(2p - l)e 2ie + e Aie 
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and as a result we see that, for instance, the quantities 



Pm = ^ \F p ^(6)\ 2 \G p ^(d)\ 2M - 2 d9 

Qm = ^ J \G p ,<p,i/)(S)\ d9 

depend only on M and p. It follows that the exit probabilities for the Hadamard walk would have 
been exactly the same had we taken U = ^= ( j ll° r ^ = -^0 1 1 ' ms * ance ' ra ther 
than the Hadamard transform. 

For arbitrary starting states, the exit probabilities may depend on <f> and V i n addition to p, 
but this change can be compensated for by considering slightly different starting states. Consider 
the starting state |0) (a \L) + (3 \R)) and transformation U p ^^. Then the exit probability for this 
walk is 



»2tt 



r M (a,(3) = *\aF PtM (0) + PG p ^(e)\ 2 \G p ^(e)\' 



But since fp,^ = f P ,o,o and e % ^ + ^ g P ,<j>,ip = g P ,o,o, we see that the exit probability is 

precisely the same as the exit probability for the walk given by unitary transformation C/ Pi o,o and 
starting state ae^'^ \0,L) + pjl-W) \0,R). 

Consequently, it suffices to study the simpler type of transformation U Pt o t o, i.e., transformations 
of the form 

\L) -^^p\L) + y^~p\R) and \R) -> ^~p \L) - ^p \R) (10) 

for p G [0, 1] to determine the properties of more general walks. Such transformations have been 
considered by Yamasaki, Kobayashi and Imai [El]]. In this case we get 



fu(z) 
9u{z) 



1 + z 2 - Vl + 2(2p - l)z 2 + z 4 



1 - z 2 - y 7 ! + 2(2p - l)z 2 + z A 



Then 



1 f 27r 

Pm = r- / \F v {9)\ 2 \G v {e)\ 2M - 2 de 

27T Jo 



1 







Letting 6 p = — — ^ 2p ^ , and T p = (—6 p ,9 p ) U (ir — p ,tt + # p ), we see (after some algebra) that 
\G V (6)\ = 1 for 8 € T p and |G ff (fl)| < 1 for 9 T p . This gives 

I / d , = coe-'(l-2p) = S m-'(2,-l) + l 



2lT r 7T 7T 2 
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This agrees with Conjecture 1 of |21] and Eq. (^Tj) below, obtained by the eigenvalue method (since 
cos (1 — 2p) = 2 sin -1 y/p). We also have 



Pm ■ ; + jz 9 > COS-^l - 2p) - ~j~ ~~ r • (12) 

This can be computed as 1/(2-71") f r \F/G\ 2 d6, but we leave this integration job to the interested 
reader. This value also follows from Eq. (|3l| ) below. 



4 Eigenfunction method 

In this section we present an eigenfunction method for computing absorption probabilities of quan- 
tum walks. We present the method for the quantum random walk introduced in Ref. [[H]]; the 
calculation for the general quantum walk corresponding to the transformation of Eq. (||) is straight- 
forward but involves slightly more complicated notation. Letting L(n, t) denote the amplitude of 
state \n,L) at time t and R(n,t) denote the amplitude of state \n, R) at time t, the dynamical 
equations are 

(L(n,t)\ = (y /p L(n + 1, t - 1) + y/l^p R(n +l,t- 1)\ 
\R(n,t)J \^T^p- L(n-l,t-l)- y/p R(n-l,t-l)J ' 

where we are ignoring boundaries for now. The Hadamard walk corresponds to the choice p = 1/2. 

Though our main interest here is in studying systems with one or two absorbing boundaries, it 
will be very useful to have at our disposal the eigenfunctions for systems with no boundaries and 
for periodic systems. Therefore, we compute them first. 



4.1 Systems with no boundaries 

If one looks for solutions of the form 



L(n,t)\ _ (A k \ i(kn-u> k t) 
R(n,t)J \B k > 



then to satisfy Eq. ( |13| ) one must have 



with 



y/p e lk y/l-p e 



ik 



^V^P e~ ik -y/p e~ ik 
The characteristic polynomial of Uk is 

X 2 - y/p [e ik - e" ifc ) A — 1, 

so its eigenvalues are 

\ k ± = y/p i sin k ± \/cos 2 k + (— 1 + 1/p) 



(13) 



(14) 
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Since Uk is unitary we may write X^± = e lu3k± with 

u>k+ = — sin -1 (y^sin k) 
Uk- = 7T — u k+ . 



(15) 



The corresponding eigenfunctions have the form (13) with 

1 



A 



k± 



<2N 



1± 



cos k 



\J\j p — sin 2 k 



B 



k± 



'2N 



1 T 



cos A; 



y/l/p — sin 2 



(16) 



We have normalized these eigenfunctions so that |Afc iCr | 2 + |I?fc i(T | 2 = 1/iV, which makes the proba- 
bilities over any iV consecutive lattice sites sum to 1. We have also chosen the arbitrary phase to 
ensure that the A k ± are real and positive. The group velocity is 



v± 



duj± 



cos k 



(17) 



V- is positive for —ir/2 < k < vr/2, and u+ 

7T/2 < k < 7T. 



y/l/p — sin 2 

= — v_ is positive for — ir < k < —ir/2 and for 



4.2 Model with periodic boundary conditions 

We now consider a system whose boundary conditions are periodic with period N. Any wavefunc- 
tion of this system can be written as a linear superposition of the eigenstates computed above: 



L(n,t) 
R(n 



:?)) - £ ^ [(&) 



For an initial condition of the form 



the Ck± satisfy 



fc£( — 7T ,7r) 

<T=± 



L(n,0) 
fl(n,0) 



ikn—iix>k a t 



(18) 



3«0 



iV 1/3 



Solving for the yields 



(19) 



In the sequel we will consider periodic systems with N large compared to the physical feature of 
interest (i.e., M <C N), with the idea of letting N — > oo. As an alternative to this procedure, one 
could also average over the continuous variable A; from the start. 
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"± ^ S JL \t > ► R 

M 

Figure 1: Diagram of system with unitary time evolution with dynamics that for all n < M are 
identical to those of a quantum walk with an absorbing wall at position n = M. 



4.3 Model with one absorbing boundary 

Now consider a system with an absorbing wall at location n = M, where for definiteness we will 
take M > 0. For our purposes it is useful to think of the wall as a boundary through which right- 
movers can be transmitted, so that the problem remains unitary and probability is conserved. We 
extend the problem so that inside the domain the dynamics are identical to the original model, 
and outside the domain right-movers move right and left-movers move left, as shown in figure |l|. 
Specifically, the model is 



n < M - 1 
n> M 



L{n,t)\ (y fp L(n + 1, t - 1) + V / T^p R(n + 1, t - 1) 

R(n,t)J Vv 7 !^ L(n - l,t - 1) - y/p R{n - l,t - 1) 

L(n,t)\ (L(n + l,t- 1 

R(n,t)J ~ \R(n-l,t-l 



together with the boundary condition 

(L{M - 1, t)\ _ ( L(M, t - 1) 

\R(M -l,t)J ~ \VT^p L(M -2,t-\)-y/p R[M — 2,t — l] 



(20) 



There are no left-movers outside the domain, and once the right-movers go through the barrier they 
are no longer converted to left-movers, so one must have L(M, t) = at all times t. The boundary 



condition Eq. (20) then implies that L(M — 1, t) = also at all times t. 

Thus, our system evolves according to Eq. (|l^) and must satisfy L(M — 1, t) = at all times t. 
The eigenstates of the periodic system discussed above clearly do not satisfy this condition. The 
only way that L(M — l,t) can vanish at all times t is for contributions from different k that have 
the same value of uj to interfere destructively (otherwise the contributions can cancel at some but 
not all times). From the dispersion relation Eq. (|i~4[), we see that u>k± = u)u—ia±, and that there 
are no other degeneracies. Thus one is led to look for eigenfunctions of the form 



/ Ck±(n,t) 
\K k ±(n,t) 



Mi 



k± 



Ak± \ ^ikn _|_ £ (A^- k) ±\ e i(TT— k)n 



e -»wk±t ; (21) 



where the A4± are normalization constants. The coefficients C,k± are fixed by requiring C{M— 1, t) = 
for all t, or 

A k± e ik(M -V + Ck±A in „ k)± e^- k K M -V = , (22) 
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Probability 
0.16 r 




Position n 



Figure 2: Probability \L(n,t)\ 2 + \R(n,t)\ 2 function of n and i of a one-dimensional 

Hadamard walk with one absorbing wall started from the initial condition L(n, 0) = and 
R(n, 0) oc exp(— n 2 /100) cos(7rx/10), which is overwhelmingly composed of values of k very close to 
7r/10. Two wavepackets propagate with group velocities ±v 0.689 (Eq. |l7|; both +v and — v are 
seen because the two bands propagate in opposite directions for a given k); one reflects from the 
absorbing boundary at position n = 100 with reflection probability P r ~ 0.184 (Eq. 



yielding (using the Ak± determined in Eq. (16)) 



I cos 2 fc cosfc \ i(7r _ 2 fc)(M-l) 

This situation is analogous to what one finds when one considers the scattering of a particle by a 
potential step, as discussed in many elementary quantum mechanics texts (see, e.g., |fT7f ; also see 
PI ). A rightmoving wave hits the step; the reflected wave is leftgoing. Figure |2[ which shows the 
time evolution of a Hadamard walk starting from an initial condition that is a superposition of a 
small band of k's, shows that an absorbing wall indeed acts in this manner. 

We know that in the — band the wavefunctions proportional to e tkn with wavevectors k in the 
range (— 7r/2,7r/2) are rightgoing, and we interpret with component at wavevector tt — k as the 
leftgoing piece generated by reflection off the boundary. The probability that a wave is reflected, 
P r (k), is just (for k G (-vr/2, vr/2)) 



„ ... l9 / / cos 2 k cos A; \ , . 

P Ak) = K ^ = ^l + -^ p - 7 = m =) . (23) 

Similar calculations for the reflection coefficient of rightmoving wavepackets in the other band and 
also for leftmoving waves that are reflected at a boundary at which R(—M + l,t) = at all times 
t yields that the probability of reflection at wavevector k in all cases is 



rw,x i> ,9 ( L cos 2 k Icosfcl \ 

P rWH& _p = M 1 + ___ 7 L') . (24) 
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These reflection coefficients agree well with the results of our numerical simulations. 

We now need to write the initial condition as a superposition of eigenfunctions. Perhaps the 
simplest way to do this is to use the method of images. Thus, we consider a system with no 
boundary; one introduces an image outside the domain that is adjusted to enforce the appropriate 
boundary condition, which here is L{M — 1, t) = for all t. We consider initial conditions in which 
the "physical" particle is at the origin, 



L(n,0) 
R(n,0) 



The form of the unnormalized wavefunction derived above suggests that it will be useful to consider 
together the pairs of wavevectors k,n — k. Thus we write 



E 



fce(-7r/2,7r/2) 
CT=± 



c h 



A k a\ ikn i /~i ( ^7r— k.a \ i(-K—k)n 

e + Un-k,cr [ D I e 

£>k,a J V-D7r-fe,cr 



where the C k)Cr are given in Eq. (|19|). We now attempt to place an image so that the condition 
L(M — 1, t) = holds at all t. We guess that the image particle should be at n = 2(M — 1) (again, 
this is suggested by the form of Eq. and try writing 



L(n,t) 
R(n,t) 



E 



fcG(-7r/2,7r/2) 



C ka I p 



Ak,a \ J,kn 



e + C^-k 



A 



Btt- 



TV—k,< 



n—k,cr j ^4(77— fc)n 



I ' \£>k,aj ' \^Tv-k,aJ 



(25) 



where once again we have used the fact that u k)CT = ^n-k,a- Again, since L(M — l,t) = for all 
times, we must have for each k S (— tt/2, tt/2): 

+ D k ,aA k ,<t& ^ ^ + ^7r— k, a-^TT— fc,<r e ^ ^ = . 



It is straightforward to verify that this equation is satisfied if we choose 

Dk,cr = 



A k 



Hence, inside the domain (n < M) the time-dependent wavefunction is 



L(n,t) 
R(n,t) 



E 



,/ ) I A k , a \ e i kn:Fka+ f^~k,a \ e i^-k)ng k ^ 



fce(-7r/2,7r/2) 



Bk,a 



B-K-k,, 



(26) 



with 



Gk,cr 



n _ J(7r-2fc)(M-i)^ Ar-fc,o- 

^k,a " ^Tr — k,<r , 

A k ,a 



•k.a 



-i(7r-2fc)(Af-l) 



C k a 



A k ^ a 

A-jr- — h_f 



(27) 
(28) 
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As expected, this wavefunction satisfies Eq. (p2[). 

The wavefunction of Eq. ( ^6|) is a superposition of plane waves. In the limit of long times, the 
only components that are in the physical domain are the leftgoing waves, which for k £ (— 7r/2, tt/2) 
are ((&,+) and (ir — k,—). Therefore, Am, the probability that the particle escapes to n — > — co 
when the absorbing wall is at M, is 

a m = Yl \^ k ,+( M )\ 2 + \Sk,-( M )\ 2 ■ 

fce(-7r/2,7r/2) 

Using Eqs. (§|), (|l9|), and (0), we see that £ fc _(M) = .F fc>+ (M), so 

A M = 2 ^ l^,+ (M)| 2 , 

fce(-7r/2,7r/2) 

with 

Converting the sum to an integral, and noting that the integrals of terms that are odd in k vanish, 
we find that 

1 

A M = - I dk C M (k) , (29) 

-tt/2 



7T 



with 

C M (k) = -. , ) ( 1/ ^+ CO f fc) + cos(vrM) cos(2MM - 1)) 

\ y/l/p-sm 2 kj V HP- 1 

+ r 2 ( 1 C ° Sfc j (1 - cos(vrM) cos 2 (fcM)) (30) 

y y/l/p - sin 2 kj 

+ 2^r cos $ [ / C ° Sk | | 1 C ° Sfc : | (cos fc - cos(vrM) cos(/c(2M - 1))) , 

V V 1/P - 1 / V VVP - s in k J 

where we have written a = le 1 ® 1 , f3 = re** r , and = — <j) r . We obtain 
Am(p) = ? 2 



\ vryTT^T l-p\ vr VA v 



r /2 dk cos L ^ COS (2fc(M-i)) 

7T 7-7T/2 Vl/p-sin 2 A: 



t/2 \J\j p — sin 2 A; 



+ r^l-^sin-^+^-l)^ 



dk =cos(2fcM) (31) 

/2 V 1/p - sin 2 fc / 



+ Ir cos $ ( ; 1 - - ( 1 H 1 2/3 sin^^/pl 

+ ' / rife cos(fc(2M-l)) 

vr y/l/p - 1 7-7T/2 V 1/p - sin 2 fc 
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M 


Ce 






1 


1 - ^ w 0.36338 


1 - 1 Ri 0.36338 

7T 


2 - ^ w 0.72676 

7T 


2 


2 - - « 0.72676 

7T 


3 - - r; 0.45352 

7T 


3 - - Ri 0.45352 

7T 


3 


4 - — w 0.816901 


13 - Ri 0.479811 


11 - Ri 0.38967 


4 


14 _ 124 ~ 0.843191 


65 - SOS » 0.489196 

■57T 


53 - » 0.372765 


5 


66 - ^ ss 0.852577 


341 x ^ 46 « 0.493304 


277 - ±|P ~ 0.367488 


oo 


| - 2 » 0.86338 


l 

2 


1 - 2 » 0.36338 

7T 



Table 1: Coefficients characterizing the probability of escape to — oo of the Hadamard walk started 
in the initial state a |0, L) + /3 |0, R) by an absorbing wall at M. Here, a = ie 1 ^ 1 , (3 = re^, and 
= — <^ r . The quantity Am, the probability of escape to — oo when the wall is located at M, is 
given by Am = i 2 C t (M) + r 2 C r (M) + ir cos $C £r (M). The probability of absorption by the wall is 
1- Am- 



If instead of Eq. (|i"3|) we use the more general transformation of Eq. (0), one finds that the only 



change in Eq. (pllj is that cos<I> — > cos($ + 20). Thus, as noted in Sec. |3.4| , the absorption 
probabilities do not depend on rj or ip, and the <p dependence can be removed by suitable adjustment 
of the initial condition. 

Eq. ( |3i"|) only applies when M > 1 , because when M = 1 the initial condition is inconsistent with 
the boundary condition unless a = 0. This complication for M = 1 is easily handled by evolving 
the walk for one time step by hand; for the initial condition a\L, 0) + /3\R, 0) the probability of 
escape to — oo when M = 1 is related to A2L, the probability that a walk starting in the state \L, 0} 
escapes to n — > —00 when the wall is at M = 2, by 



Ai = lay^ + ZV 1 "Pi A 2L • 

For the special case of the Hadamard walk (p = 1/2, <j) = 0), we find k\ = (1 + 2^r cos $)(! — 2/-7r). 
Values for other finite values of M are easily obtained by evaluation of the integrals in Eq. ; a 
few values for the Hadamard walk are listed in Table |l]. 

In the limit M — > 00 all the M-dependent terms vanish because of the oscillating integrands. 
The coefficient of the r 2 term in the limit M — ► 00 has been correctly conjectured by Yamasaki 
et al. on the basis of numerical results. In the limit M — * 00, the probability that a particle 
undergoing a Hadamard walk escapes to n — > —00 is 

A ~ = /, (i-i) + T +fr »*( 1 -|) ■ <32> 

As expected, this agrees with Eq. (H). Integration by parts of the integrals in Eq. ( |3i~| ) may be used 
to find the behavior at large but finite M (c.f., Eqs. (|5|- 



4.4 Model with two absorbing boundaries 

Now we discuss the situation in which there are two absorbing walls, so the domain of the quantum 
walk is finite. We restrict consideration to the Hadamard walk (p = 1/2, 77 = ip = (f> = 0), but no 
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essential complications are introduced when rj, <j>, ip, or p take on different values. We place the 
walls at +Mr and —Ml and again start the walker out at position n = 0. 

By exploiting the results obtained above for the single wall case, it is rather simple to compute 
the absorption probabilities when Ml is very large. This can be seen by noting that because when 
the walls are very far apart, the evolution can be viewed as a succession of reflections off the two 
walls followed by a single transmission. From Eq. (^4|), we know that the reflection probability is 
P r = (\/T+ cos 2 k — | cos/c|) 2 . Let us denote by l~LM R {k) the probability of transmission through 
(or, in other words, absorption by) the left boundary of the component at wavevector k when the 
right wall is at Mr. Because the wavepacket can be absorbed at —Ml the first time it hits that 
boundary, or it can be reflected at both the left and right boundaries and then absorbed at the left 
boundary, etc., we have 

T LMR (k) = C MR {k)[{l-P T ) + P r P r (l-P r )+P r P r P r P r (l-P r ) + ...} (33) 
= ^{k)^- (34) 
1 „ , , x / cos k 



where £M R (k) is defined in Eq. (3C). 

When Mr > 1, the transmission through the left wall when the right wall is located at Mr, 
Tlm r , is thus 

1 /-t/2 

T L m r = - / dkT LMR (k) , (35) 

J-tt/2 

with 

TLM R (k) = -I 2 - (2 + cos(2A;) + cos(ttM) cos(2A:(M - 1))) (36) 

2i L ~\~ COS rh 

r 2 1 

+ 5t(! " cos(ttM) cos(2A:M)) (37) 

2 1 + cos z k 

+ Ir cos $ cos k 1 (cos k - cos(ttM) cos(k(2M - 1))) . 

1 + cos^ k 

Once again the Mr = 1 case is done by evolving the system by hand for one time step. The results 
for some finite values of Mr are listed in Table 2. The result for Mr = 1 was obtained previously 
in Ref. §. 

As Mr — > oo, the absorption by the left wall, Tloo is 

The minimum value of Tloo is when $ = 7r, t = sin(7r/8), where T Loo = 1 - 1/^2 xi 0.292893, and 
the maximum value of Tloo occurs when $ = 0, I = sin(3vr/8), where T Loo = l/y/2 w 0.707107. 

When the two walls are close together, the calculations are straightforward in principle but 
tedious in practice. The overall strategy of the calculation in this case is exactly the same as in all 
the calculations above — instead of considering a problem with absorbing barriers, one considers a 
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M 


Di 






1 


1 7= ps 0.292893 

V2 


1 7= ps 0.292893 


2 - V2 ps 0.585786 


2 


2 — v2 ~ 0.585786 


6 — 4y 2 ps 0. 343146 


z~* a / c\ r\ 0(01 a /"> 

6 — 4y2 ps 0. 343146 


Q 

o 


f — -y= ~ U.DooUoy 


00 — -y= PS U.OOl/DO 


on 01 /o ^ n Qm £1 £ 
oU — zlyz ps U.oUioio 


4 


36 - 25\/2 ps 0.64461 


204 - 144V2 ps 0.353247 


170 - 120^/2 ps 0.294373 


5 


205 - J| ps 0.64614 


1189 - ^§ ps 0.353501 


986 - 697\/2 ps 0.293147 


oo 


1 - 5J5 « 0.646447 


A- « 0.353553 

2v 2 


1 - 4= ~ 0.292893 



Table 2: Coefficients characterizing the probability of escape through the left barrier to -co of the 
quantum walk started in the initial state a |0, L) + (3 |0, R) with absorbing walls at M and at —Ml, 
when Ml — > 00 Here, a = £e lr ^ 1 , (3 = re l< ^ r , and = — <j) r . The quantity 7l(M), the probability 
of escape to -00, is given by T L (M) = £ 2 D e (M) + r 2 D r (M) + trcas®D lr (M). The probability of 
absorption by the right wall at M is 1 — Tl(M). 



unitary problem whose time evolution is identical within the physical domain. Here we embed the 
walk in a large periodic domain; the model is 



M L + 1< n < M R - 1 



L(n,t)\ 1 (L(n + l,t-l)+R(n + l,t-l] 

R(n,t)J ^/2\L(n-l,t-l)-R(n-l,t-l] 



^ , ^ ( 'L(n,t)\ f L(n + l,t - 1)\ 

n>M R ^n<-M L : ^J } j = [ R \ n _ ,\ f _ ,]) (38) 

(L{M R -l,t)\ ( L(M R ,t-l) \ 

\R(M R -l,t)J (L(M R - 2, t - 1) - R(M R -2,t-l))) 

L(-M L + 1, t)\ ( ^ (L(-M L + 2, t - 1) + i?(-M L + 2, t - 

i?(-Mi-l,^ ^ R(-M L ,t-l) J ' 

We wish to find the energy eigenstates that satisfy 



L(n,t)\ = (Luin) 
R(n,t)J \Ru{n) 



and are consistent with the boundary conditions L(M R + N) = L{—M£), R{M R + N) = R{—Ml)- 
The eigenvalues and eigenfunctions of the model are specified by these conditions and can be 
found similarly to those of systems with a finite square well potential |jl3[| . In the "outer" region, 
the eigenfunctions take the form 



-iuit 



where the periodic boundary condition is enforced by requiring e tu) ( M L+N) _ e -iuM R .^ ^ e u anc j 
v will be fixed by matching to the inner region. Since we wish to take N —* 00, there will be a 
continuum of values of u. Finding the wavefunction for all Ml and M R is rather involved; for 
every u one must match the inner and outer wavefunctions, keeping in mind that in the inner 
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region some of the wavefunctions are of the form e Kn with k real (this is because the quantum 
walk has a bandgap). When the inner region is small, the eigenfunctions can be found by brute 
force matching to the outer region. For Ml = 1, Mr = 1, we have L w (0) = and i? w (0) = v u . 
The wavefunction must satisfy the initial condition, so = a and v u = [5. Thus, for this case, 
as expected, the probability of absorption at the right and barriers are \(3\ 2 and |a| 2 , respectively. 
When Ml = = 2, the coefficients Aj(n) in the "inner" region satisfy 

Mi) 

L u (0) 

These equations can be solved to yield L w {ti) and -R^(l) in terms of w w and v u , which are in turn 
fixed by enforcing the initial condition. However, the algebra is complicated and will not be given 
here. 



—6-^(0) -Km 



(39) 



4.5 Arithmetic properties of exit probabilities 

In this subsection we discuss what can be said about the exact exit probabilities for the Hadamard 
walk. 

Examination of Table 1 above suggests that the coefficients Q(iVf), C r (M), and C^ r {M) are all 
rational numbers plus rational multiples of 1/ir. This can be proved as follows. Taking p = 1/2 in 
Eq. @ we get for M > 1 

Cl (M) = ( 3 — 2 -)-t±}l f" 2 co S kc OS (2(M- m ^ 



c, m -i + t$L r' 2 c ° s < 2Mfc ) iK 



2 TT J TT J-. n /2 V 1 + COS 2 k 

2 ' TT J-n/2 VT+ COS 2 k 

and 

c, A M) = l -l + h-i) M cos2fccos((2M '" t) *. 

TT J -it/2 V 1 + COS 2 k 

In the notation of section 3, these are 1 ~Pmi 1 — Qm > an d — %(pq)mi respectively. 

Using the identity e imk = (cos k + isink) m and the binomial theorem, we can express cos(mk) 
as a polynomial in cos k. Furthermore, this polynomial has only odd powers of cos k when m is 
odd, and only even powers for m even. Making the substitution u = cos k, we see (since we started 
with even integrands) that each of the integrals is a linear combination, with rational coefficients, 
of expressions of the form 

v.* du _V^ r((j + 1)/4) 



VT^? 4 r((j + 3)/4) 

with j odd. For such j, one of {(j + l)/4, (j + 3)/4} is an integer and the other is a half integer. 
Thus one of the T values is an integer and the other is a rational multiple of y/n, giving the desired 
form. 
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A similar argument can be used to prove that the coefficients Di, D r , and De r of Table 2 are 
rational numbers plus rational multiples of y2. For this one needs the integral 



u 2 i 



du. 



lo (1 +u 2 )Vl -u 2 
If we substitute u 2 = t and use Theorem 2.2.2 of we see it equals 

r( J + i/2) 



2 r(j + 1) 



2F1 



1, j + 1/2 



where 2F1 is the Gaussian hyper geometric function. It then follows from (15.3.5) and (15.4.11) of 
that the coefficients have the desired form. 

These observations allow us to generalize a result of ||. Suppose we start both walks at the 
same distance from the right barrier, in the state a = 1, (3 = 0. Then the escape probability for the 
1-barrier walk is never equal to the limit (as the left barrier recedes to —00) of the corresponding 
probability for the 2-barrier walk. This is so because ir is not an algebraic number. The physical 
interpretation of this result is that no matter how far away the second wall is, it reflects probability 
back to the first wall. 

The precise arithmetic properties of these numbers are interesting to contemplate. We can show 
by a combinatorial argument, for example, that when Ce(M) = cmi + c-miI^i the denominator of 
cm,2 divides the product of the odd numbers in {1, . . . , M — 1}. It also appears that cm,i is always 
an integer, but we have not yet proved this. 

4.6 Time dependence of absorption 

The eigenfunction method can be used to calculate the time evolution of the absorption by the 
boundaries. We consider here the situation when the walker starts off far from any boundary 
(M large). We also specialize to p = 1/2; once again, generalization to other transformations is 
straightforward. 

When there are two walls, the approach to the asymptotic behavior is quite slow. This is 
because the reflection coefficient tends to unity as k — ► 7r/2, and the group velocity v g also vanishes 
as k — * 7r/2. Therefore, it takes a very long time for all the probability to be absorbed at the 
walls — the fast-moving wavepackets are absorbed quickly, but as the evolution proceeds one is left 
with components that move slower and slower and are absorbed less and less. The fraction of the 
probability at wavevector k that is absorbed per unit time is 



[fraction of probability absorbed per collision] 



2 

1 — ( \/T+ cos 2 k — cos k . 



[time between collisions] 2M/v g 
Thus, we estimate that the probability at wavevector k obeys the differential equation 

dV k 



dt 



-lkP k {t) , (40) 



with 



1 — (\/T+ cos 2 k - cos k) 2 , . 

lk ~ 2MVlTms r k/ cosk ' 
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Figure 3: The time dependence of V(t), the probability that the particle has not been absorbed, is 
different in systems with one and two absorbing boundaries. With one boundary, V(t) oc V{oo) + 
constant ft 2 (Eq. (p 



])), while with two boundaries, V{t) oc 1/yi (Eq. (42)). The faster approach 
to the asymptotic value for systems with one boundary is because there are no repeated reflections, 
while the constant probability at short times for both systems is because no probability has yet 
reached the boundary at position M. Simulations and the analytic forms are all obtained using the 
initial conditions i = and r = 1. The probability at t = oo for one boundary is given by Eq. (j31~l), 
and for two boundaries it is zero. 



so that Vk(t) = "Pfc(0)e _7fef . The total probability that the particle is inside the region at time t, 
'P(i), is the integral over wavevectors k of the probabilities Vkit). At long times the probability 
is dominated by /c's near ±7r/2 (since the decay rate tends to zero there); expanding for fc's near 
±7r/2 yields 

V(t) = - dk e -(^m 2 t/M 

n Jo 

Now changing variables to s = (k — 7r/2)\/£ and letting the limit of integration \\ft — > oo gives 

V{t) = —r if ds e-VM = 

When there is one wall that is far from the origin, the approach of the absorption probability to 
the asymptotic value is faster because there is no possibility of repeated reflections. The behavior 
can be characterized using Eq. ( |25| ) by calculating the total left-going probability from the physical 
particle, 

r' 2 dk (ic fc ,+i 2 +ia_ fc ,„i 2 ), 

J-tt/2 

the fraction of right-going probability from the physical particle that has yet to reach the wall 
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position n = M, 

»7r/2 



r , dk(\c k , + \ 2 +\c^. 

icos-M -. m ) 



fc, 



and the fraction of left-going probability from the image particle that has already passed the wall 
position n = M, 

_1/ M/t \ 

/cos ( ) 

as a function of time. The C"s and -D's represent probability amplitudes of plane waves emanating 
from the physical and image particle, respectively, as indicated in Eq. (25), and the limits of 
integration take into account that at time T, only the fc's for which v g T < M have yet to reach the 
wall. This gives the asymptotic time dependence for the probability as 

M 2 M 3 

p( t)= p(oo) + -^+0(- r ), (43) 



where V(oo) = Am{p = 1/2) is given in Eq. (pip. Figure || demonstrates that the asymptotic forms 
(|42]) and (p3| ) provide an excellent description of the time dependence for times greater than a few 
times the system size. 



5 Quantum walks in general dimensionality 
5.1 Introduction 

A quantum walk on a d-dimensional lattice n G Z d may be defined in several ways. We shall adopt 
the following. On each site there are "internal degrees of freedom" labeled by an integer s £ [1, 2d], 
where s replaces the L, R notation of the previous section. The time-dependent wavefunction 
is written as ^(n,s,t). The index s\ corresponds to a left move in the first dimension: n = 
(ni,n2~-nd) — > n = ( n i — l ; ^2"- n d) 5 £2 corresponds to a right move in the first dimension: n = 
(ni,7i2...7i(/) — ► n = {n\ + l,n2...n^), S3 corresponds to a left move in the second dimension: 
n = (ni,n2~-Tid) ~ > n = ( n ij n 2 — ^--- n d), an d so on. Each time step consists of the "coin-toss" step 

2d 

V(n,s,t) -» C ss rV(n,s',t) = y'(n,s,t), 

s'=l 

where C ss / is a real 2d x 2d orthogonal matrix, followed by the walk step 

V(n, s,t) ->• ^ + ^ s -5' *) = s > * + 1 )- 

Here 5 £ Z d is a nearest neighbor separation: |<5| = 1. 6 takes on 2d possible values, sg is the 
index corresponding to a move in the ^-direction. This sequence of two successive transformations 
defines a unitary operator U^{n,s,t) = s,i + 1). 

We shall assume that there is a separate coin toss for each dimension. This makes C ss i into a 
block-diagonal matrix consisting of d real 2x2 orthogonal matrices along the diagonal and zeroes 
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elsewhere. A more general definition allows all the elements of C ss i to be nonzero. Yet a third 
definition is that adopted by Mackay et al. JTl]]. These authors introduce d qubits at each site, 
resulting in a 2 d x 2 d coin-toss matrix. The second and third models can also be treated by the 
eigenfunction method, but the results are more cumbersome. 

In this section, we give the formal eigenfunction expansion solution for the asymptotic behavior 
of our ci-dimensional quantum walk (first model) and the survival probability in this model when 
ad — 1-dimensional absorbing wall is present. 



5.2 No boundaries 



The problem of finding ^(n,s,t) = J7*\I/(rt, s,0) with given ^(n, s,0) in the absence of absorbing 
walls may be solved by an eigenfunction expansion. We make the Ansatz 



*(ra,s,i) = t/jq{s)e 
where ipq(s) are the eigenfunctions of U~ s ': 



i(q'ft—uifft) 



f V>«r(s)> 



(44) 



and the matrix Ui s is given by 



( Cue** 




V - 



C 22 e- iqi 








Cue-** 






Cue** 
Cue~ iq2 



J 



The momentum variable q = (gi, q 2 , ...qd) lies in the hypercubic region \qi\ < n, i = l...d. (It 
is convenient to envision the hypercubic lattice to be extended periodically in all dimensions, and 



then taking the period to infinity.) There are in fact 2d solutions to the eigenvalue equations 44 at 
a fixed value of q. Let us index the eigenvalues u>^ and the eigenfunctions ^2 by a = 1...2d. This 
band index a replaces the (+, — ) labels of the previous section. 

The eigenvalues of the operator U have important degeneracies. 

Let us first consider the 2x2 submatrix u\: 



it. 



(Qi) 



Cue 
C 21 e 



C 12 e 
C 22 e 



"/i 



-iqi 



This is unitary and satisfies detwi = ±1. A determinant of —1 for a submatrix characterizes 
the Hadamard transformation, but a general 2x2 coin-toss operator could also have a positive 
determinant. One can show that u\ satisfies the identity Rui^q^R -1 = det(ui)ui(—qi), where R 
is the matrix 



R 



-1 

1 



Thus u\(qi) is unitarily equivalent to u±(— qi) up to a sign. If we denote the reflection (qi,q 2 , ■■■qd) 



(-qi,q 2 , ...qd) hyq^ q R we have u/- 



1R 



if det Ui (q) = 1 and 



-9l>92,-9d) 



if det u\{q) 
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For a general coin toss matrix C ss i no such simple description of the degeneracies exists. 
A quantum walk that starts at the point n = (0, 0, ...0) with arbitrary internal state satisfies 

v(n, s,t = o) = %*o(«) = / E a ^ ^ em °^ ( 45 ) 

so that the coefficients a a are determined by the initial conditions. We shall employ the shorthand 
notation 

and note that 

/ J*** = W 

Thus, given the initial condition, we find the a a (q) by inverting a set of linear equations at each 
value of q: 

a 

The full solution is then 

a J 

To understand the asymptotic properties of this wavefunction at long times, consider a point that 
moves at constant velocity n = ct away from the origin. After this substitution, the behavior of 
the integral as t — > oo may be obtained by the stationary phase method. The neighborhood(s) of 
the point or points in g-space (n) where — c = dominate the integral. The result is 

a ^ 

d/2 

x £5> Q (#)^(s)exp [(iqf ■ c- iu%)t] |Jf |~ 1/2 t^ 2 + ©(t"^ 1 ), 

a i 

where J™ is the Jacobian of the function u;2 at the point qf . If there are no solutions to the 
zero-gradient condition, then the asymptotic behavior of \^(n,s,t)\ is generically determined by 
contributions from the boundary of the region of integration and one finds \^(n, s,t)\ = 0(t~ d ). 

From this expression, we deduce that the probability |*(n, s, t)\ 2 spreads linearly with time. The 
falloff of the wavefunction as t~ d l 2 implies the following physical picture. The initial wavepacket 
is a superposition of waves with various group velocities = V^lo^. Each such component moves 
according to the ballistic equation n = v^t. Because of the limited range of q, there is a maximum 
group velocity in each spatial direction. This maximum velocity defines the wavefront in that 
direction. All components move at constant speed and the overall probability is normalized: 

^l^s^)! 2 ^^ x est. ~ r d t d -► l. 

its \n\<t 
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5.3 (d — 1) -dimensional absorbing wall 



We now generalize the absorption problem to a (d—1) -dimensional wall located at n = (M, 0, 0, 0) 
with M > 0. We shall treat only the case where detiii = +1, as the other case has been treated in 
detail in one dimension. 

In fact, many of the results from See]!] generalize immediately. We again extend the problem 
to the full space, stipulating that there is no motion to the left in the region n\ > M. Then it is 
sufficient to solve for the wavefunctions ^f(n,s,t) = ^f(n\,n2, ■■■rid, s,t) in the region n\ < M that 
satisfy ^(M — l,ri2, —Ud, l,t) = 0, which means we have solutions of the form 



2d . 

VP(n, s,t) = J2 K(£ M) m{8) + A a (q R , M) m(s) 

i Jqi>0 L 



e 9 Dq. 



a=l 

Enforcing the boundary condition leads to the relation 

2d 2d 



a=l a=l 



between an expansion coefficient and its reflected counterpart. As seen already in Sec.||, this relation 
is consistent with the initial conditions, if it is kept in mind that the wavefunction is only needed 
in the half-space ri\ < M. The initial condition is 

2d . 

9{ft, s,t = 0) = S H ^ (s) = W A a (q, M) ^(s) e^Dq. 

Using the method of images or otherwise, we invert this expression to determine the A a (q, M) and 
the final solution for n\ < M is 

2d . 

To obtain the survival probability Am, we note that at long times only leftmoving waves (vi(q) = 
duj^/dqi < 0) will be present in the physical domain: 

2d 2d 



A M = V V f \A a {q,M)^(s)\ 2 Dq. 

' „_i Jvf(q)<0 



a=l s=l 



This expression manifestly satisfies < Ajy < 1, except for special choices of the initial condi- 
tion. This is in sharp distinction to the classical random walk, for which the survival probability 
vanishes for all d and M. 
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